An effective Hamiltonian for an extended Kondo lattice model and a possible origin of 

charge ordering in half-doped manganites 
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An effective Hamiltonian is derived in the case of the strong Hund coupling and on-site Coulomb 
interaction by means of a projective perturbation approach. A physical mechanism for charge 
ordering in half-doped manganites (Ro.sXo.sMnOs) is proposed. The virtual process of electron 
hopping results in antiferromagnetic superexchange and a repulsive interaction, which may drive 
electrons to form a Wigner lattice. The phase diagram of the ground state of the model is presented 
at half doping. In the case of formation of Wigner lattice, we prove that spins of electrons are 
ahgned ferromagnetically as well as that the localized spin background is antiferromagnetic. The 
influence of the on-site Coulomb interaction is also discussed. 



PACS numbers; 71.45.Lr, 75.30. Fv 
I. INTRODUCTION 

The family of doped manganites, Ri-ajX^jMnOa 
(where R = La, Pr, Nd; X=Sr, Ca, Ba, Pb), has renewed 
both experimental and theoretical interests due to the 
colossal magnetoresistance and its potential technologi- 
cal application to magnetic storage devices. Apart from 
their unusual magnetic transport properties, experimen- 
tal observations of a series of charge, magnetic and orbital 
ordering states in a wide range of dopant also stimulate 
extensive theoretical curiosities. Early theoretical studies 
of manganites concentrated their effort on the existence 
of metallic ferromagnetiam. From the so-called "Double 
Exchange" (DE) modeljil in which the mobility of itiner- 
ant electrons forces the localized spins to align ferromag- 
netically, one can understand qualitatively the relation of 
transport and magnetism. However the rich experimen- 
tal phase diagrams are far beyond the DE model. For 
example, according to the DE model, itinerant electrons 
have the lowest kinetic energy in a tight binding model, 
and should be driven to form a more stable ferromagnetic 
phase when the system is half doped, i.e., x=0.5. On 
the contrary, it is insulating rather than metallic ferro- 
magnetic at a low temperature as expected theoretically. 
Furthermore, a charge ordered state was observed, which 
is characterized by an alternating Mn'^+ and Mn'*+ ions 
arrangement in the real spaceB Usually when the repul- 
sive interaction between charge carriers dominates over 
the kinetic energy the charge carriers are driven to form 
a Wigner lattice. It has been shown experimentally that 
the charge ordering is sensitive to an applied magnetic 
field at low temperatures: resistance of a sample may de- 
crease in several order of magnitude and the charge order- 
ing disappears at a low temperaturejj which implies that 
the repulsive interaction should have a close relation to 
the spin background. Although there have been extensive 
theoretical efforts on anomalous magnetic propertiesjj a 
comprehensive understanding on the physical origin of 
ordered states and their relations to the transport prop- 
erties are still awaited. 

To explore electronic origin of these phenomena, we 



try to establish a more unified picture to understand the 
physics starting from an electronic model, which has been 
used to investigate the magnetic properties of the system 
extensively. We derive an effective Hamiltonian in the 
case of the strong on-site Coulomb interaction and Hund 
coupling by means of a projective perturbation approach. 
It is found that the virtual process of electron hopping 
produces an antiferromagnetic superexchange coupling 
between localized spins and a repulsive interaction be- 
tween itinerant electrons. The antiferromagnetic correla- 
tion will enhance the repulsive interaction and suppress 
the mobility of electrons. In the half-doped case, i.e., 
X — 0.5, relatively strong repulsion will drive electrons to 
form a Wigner lattice. In the case of the Wigner lattice, 
we prove that the electrons are fully saturated while the 
localized spins form an antiferromagnetic background. 
Strictly speaking, the ground state possesses both anti- 
and ferromagnetic, i.e., ferrimagnetic long-range orders. 



II. EFFECTIVE HAMILTONIAN 

The electronic modeLfor doped manganites studied in 
this paper is defined asEl 

- Jh ^ Sj • Sic + Jaf ^ Si • Sj. (1) 

i (ij) 

where cj ^ and Ci^a are the creation and annihilation op- 
erators for Cg electron at site i with spin a (=Tji)i re- 
spectively, (ij) runs over all nearest neighbor pairs of 

lattice sites. Sic = J2a a' '^'^'^''^1 a'^i,'^' /'^ ^r*5 the 

Pauli matrices. Si is the spin operator of three t2g elec- 
trons with the maximal value 3/2. Jh > is the Hund 
coupling between the Cg and t2g electrons. The antifer- 
romagnetic coupling originates from the virtual process 
of superexchange of t2g electrons. In reality, the Cg or- 
bital is doubly degenerated. For the sake of simplicity, we 
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only consider one orbital per site, which amounts to as- 
suming a static Jahn- Teller distortion and strong on-site 
interactions (relative to kinetic energy). 



FIG. 1. Tow virtual processes of electron hopping in the 
restricted Hilbert space which favors antiferromagnetic cor- 
relation between neighboring sites. On the left side are the 
initial states, and on the right side are the mediate states. The 
process (a) leads to an effective attraction between electron 
and hole, and the process (b) leads to an effective attraction 
between electrons. 

Usually the Hund coupling in the doped manganites 
is very strong, i.e., JhS :s> t. Large JhS suggests that 
most electrons form spin S" -f 1/2 states with the local- 
ized spins on the same sites, which makes it appropriate 
to utilize the projective perturbation technique to investi- 
gate the low-energy physics of the Hamiltonian (|l|) . The 
effect of finite and large JhS can be regarded as the 
perturbation correction to the case of infinite Jjj, which 
is described by a quantum double exchange modeLQ Up 
to the second-order perturbation correction, there are 
two types of the virtual processes which contribute to 
the low energy physics (See in Fig. 1): (a). An elec- 
tron hops from one site to one of the nearest neighbor 
empty site to form a spin S* — 1/2 state and then hops 
backward. The intermediate state has a higher energy 
AE'a = J// (5 -I- 1/2) than the initial state, (b). One elec- 
tron hops from one site to one of the singly occupied sites 
and then backward. The intermediate state has a higher 
energy Ai?b = JhS + U than that of the initial sts 
Hence, by using a projective perturbation approach| 
the effective Hamiltonian is written asQ 
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Si is a spin operator with spin S + 1/2, and a combina- 
tion of spin of electron and localized spin on the same 
site. Pih and P^ are the projection operators for empty 
site and single occupancy of spin S + 1/2. The hp 
term in Eq. (0) is the quantum double exchange modelJ 
It enhances ferromagnetic correlation, and may be sup- 
pressed if the antiferromagnetic exchange coupling of lo- 
calized spin is very strong. The second, third and fourth 
terms prefer antiferromagnetism to ferromagnetism. The 
third term describes an attractive particle-hole interac- 
tion since the value of the operator before PuiPjs 
ways non-positive. In another words, an repulsive in- 
teraction between electrons in the restricted space arises 
when the spin background deviates from a saturated fer- 
romagnetic case. 

To simplify the problem, we take the large spin ap- 
proximation, and keep JhS = jh and JafS^ — jaf - The 
spin operator is parameterized in polar angle 6 and <j). In 
the approximation, the Hamiltonian is further reduced 
to 
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where 



c = cos - cos ^ + sin ^ sin rig-^t^'"^.) 
^ 2 2 2 2 



cos Qij — cos 9i cos 9j + sin 9i sin 9j cos( 
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Physically, a is an electronic operator which is fully po- 
larized along the localized spin on the same site. \cij \ = 
cos(0y /2) and approaches to zero when Qij —^ tt. If we 
neglect the Berry phase in Cy , the first term gets back 
to the classical DE model. Now it is clear that the fer- 
romagnetism is always predominant in the ground state 
if other terms in the effective Hamiltonian (Eq. (||)) are 
neglected. The sign of the interaction 
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]h 



2 2jhU + 3h 



(4a) 



is determined by the ratio jh/U. If U is less than jh, 
the interaction is attractive, but if U is greater than jh, 
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the interaction is repulsive. The attractive or repulsive 
interaction will lead to different physics. Hence U = jh 
is a quantum critical point. The influence of the on- 
site Coulomb interaction will change qualitatively (not 
just quantitatively) the physics of the doped mangan- 
ites, which is usually ignored. In the case of small U, the 
attractive interaction will drive electrons to accumulate 
together to form an electron-rich regime, i.e., the phase 
separation may occiu; when the spin background becomes 
antifpromagnetismJa Monte Carlo simulation by Dagotto 
et al£3 shows that the phase separation occurs in the case 
of [/ = 0. However, the phenomenon was not observed in 
the case of hfge U. The phase diagram of ?7 = is also 
seen in Ref.tl. From our analysis, the attractive interac- 
tion originates from the virtual process (b). Due to the 
double occupancy in the intermediate state, an extra en- 
ergy U costs in the process. When U is sufficiently large, 
the process (b) will be suppressed and the process (a) 
becomes predominant. The net interaction between elec- 
trons is repulsive. Therefore the phase separation may 
occur only if [/ < jf, . 

III. ORIGIN OF WIGNER LATTICE 

We are now in the position to discuss the instability to 
the Wigner lattice. In the doped manganites, the on-site 
Coulomb interaction is muchj^ronger than the Hund's 
rule coupling, i.e., U ^ JnS.tB In the case, the process 
(b) in Fig.l needs a much higher energy to be excited 
than the process (a) does. The process (a) dominates 
over the process (b). The effective interaction is repul- 
sive. Hence we shall focus on the case of strong corre- 
lation (i.e., U » JrS). To simplify the problem, we 
take U — > +00 and neglect the term containing U in 
Eq.(||). A finite and large U will produce minor quan- 
titative (not qualitative) changes of the physics we shall 
discuss. The ratio of the repulsion to the hopping term 
r = (i/j>i) sin^ ^/ cos ^ depends on not only t/j^, 
which is usually very small, but also the angle of two 
spins, r = if Qij = 0, and -|-oo if &ij — n. In other 
words, the ratio could become divergent in the antiferro- 
magnetic spin background (Gy — tt) even though t/jh is 
very small. Relatively large ratio will make a state with a 
uniform density of electrons unstable. To understand the 
physical origin for the Wigner lattice at x = 0.5, we first 
see what happens in the antiferromagnetic background. 
When all Qij — > tt, the average energy per bond is — 2j/i if 
the two sites are empty or occupied, and —{2jaf + / jh) 
if one site is empty and another one is occupied. The 
later has a lower energy. At x = 1/2, {a\ai) — 1/2. The 
average energy per bond is —{2jaf + /2jh) for a state 
with a uniform density of electrons. If the electrons form 
a Wigner lattice, i.e., ((aja^ - l/2)(a]aj - 1/2)) = -1/4, 
the average energy per bond is — (2ja/ + t^/jh), which 
is lower than that of the state with a uniform density. 
Therefore in the antiferromagnetic background a uniform 



density state is not stable against the Wigner lattice even 
for a small t/jh- The same conclusion can be reached 
by means of the random phase approximation. On the 
other hand, the formation of the Wigner lattice will also 
enhance the antiferromagnetic exchange coupling from 

-jaf to -~{jaf +t^/2jh). 

The phase diagram of the ground state is determined 
by the mean field approach. Several of the features are 
determined in several limits: for example, the ground 
state is ferromagnetic at t/jh = and jaf = 0. Due to the 
instability to the Wigner lattice or charge density wave 
for finite t/jh and jaf we take (a-a^ — 1/2) ~ Ae**^ ""' 
where Q = (tt, tt, • • •) and (• • •) is the ground state aver- 
age. We aJap take {aj) = cos(9/2) and (sin^(6ij/2)) = 
sin^(0/2).E3 The free energy per bond is 

^(A, e) = - / (|).^^^Wcos4+4|sin^|A^ 

1 i\ • 2 © . 2 a2 

- (jaf + T-) sm2 - + - sm^ -A2 (5) 
4j/i 2 Jh 2 

where e(fc) — — ^(X^q^i coska)/d and d is the number of 
dimension. The integration runs over the reduced Bril- 
louin zone. The phase diagram (Fig. 2) is obtained by 
minimizing the energy £^(A, O). A and & are the order 
parameters for charge and magnetic orderings, respec- 
tively. A = and = represents a full ferromagnetic 
(FM) phase, A = and O 7^ represents a canted ferro- 
magnetic (CF) phase, A = 1/2 and O = tt represents the 
Wigner lattice (WL), and A < 1/2 and O / represents 
a mixture of charge and spin density waves. A full fer- 
romagnetic phase diagram appears at smaller t/jh and 
jaf, which indicates that the double exchange ferromag- 
netism is predominant. The Wigner lattice appears at a 
larger t/jh and jaf - The antiferromagnetic coupling orig- 
inating from the virtual precess (a) and superexchange 
coupling of the localized spins can suppress the double 
exchange ferromagnetism completely. A canted ferro- 
magnetic phase is between the two phases. At jaf — 0, 
the transition from ferromagnetism to the Wigner lat- 
tice occurs at t^/jh = 2j dke{k) / (2%)'^ which equals to 
0.63662i for d = 1, 0.405282t for d = 2, and 0.336126i for 
d = 3. When the effective potential energy t/jh becomes 
to dominate over the kinetic energy, the ferromagnetic 
phase is unstable against the Wigner lattice. For a finite 
jaf, a smaller t/jh is required to form a Wigner lattice. 
However t/jh must be nonzero, even for a large jaf. In 
the double exchange model, i.e. j/i — > -l-oo, we do not 
expect that the Wigner lattice could appear at low tem- 
peratures at a; = 1/2 unless a strong long-range Coulomb 
interaction is introduced. 
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FIG. 2. The phase diagram of the ferromagnetic Kondo 
lattice model on a cubic lattice {d — 3) at x = 1/2. 



IV. FERRIMAGNETISM AND WIGNER 
LATTICE 

We go back to Eq. (||) to discuss the magnetic proper- 
ties of the ground state (or at zero temperature) in the 
case that the Wigner lattice is formed at a; = 1/2 (1 — a; 
is the density of electrons). The charge ordering in the 
manganite is an alternating Mn'^+ and Mn'*+ arrange- 
ment rather than a charge density modulation, which 
means rii = 1 or 0. A d-dimensional hyper-cubic lattice 
can be decomposed onto two sublattice A and B. In the 
charge ordering state, suppose that all electrons occupy 
the sublattice A, then 

+ _ J 1, if i G S and j e A; , . 

^rh-Tj, - ]^Q^ otherwise. 

The first term in Eq.(^ must be suppressed completely 
as the Wigner lattice is a static real space pattern, i.e., 
the hopping processes are forbidden. In the case, the 
Hamiltonian is reduced to 

where J'^p = JafS'^ + 2St'^/JH{2S + 1)^ and the sum- 
mation runs over the nearest neighbor pairs. This is an 
antiferromagnetic Heisenberg model. The spin on the 
sublattice A is S* -I- 1/2 as the electrons on the sites 
form spin S* + 1/2 state with the localized spins, and 
the spin on thCj-sublattice B is S. According to Lieb- 
Mattis theorem,liJ the ground state of Eq.(0) is unique 
apart from spin SU(2) {2Stot + l)-fold degeneracy. The 
total spin of the ground state 5*40* is equal to the differ- 
ence of the maximal total spins of two sublattices. In the 
case, 

Stot = ^ (8) 



which is also the maximal total spin of electrons (iVg is 
the number of electrons). It seems to be that all electrons 
are saturated fully while the localized spins form a spin 
singlet state. Furthermore, it is shown rigorously that 
the ground state possesses antiferromagnetic long-raage 
order as well as ferromagnetic one for any dimensionJlj 

V. DISCUSSION AND SUMMARY 

We wish to point out that, in the case that the 
Wigner lattice is formed, the magnetic structure estab- 
lished here is unlikely-4p-.full agreement with all exper- 
imental observations .tSllZl The model discuss here is a 
simplified theoretical model which has neglected some 
effects, such as the orbital degeneracy of Cg electrons, 
strong John- Teller effect and lattice distortion. Method- 
ologically, we apply the projective perturbation approach 
to deal with the model. The strong electron-electron cor- 
relations has been successfully taken into account by the 
projection process. The perturbation process tells us that 
the effective Hamiltonian should be valid at small t/jh, 
which requires a strong Hund coupling comparing with 
the hopping integral t. In practice, the parameters of the 
model for doped manganites are roughly estimated as. 
U « 5.5eV, Jh ~ 0.76eV, t « O.AleV, Jaf « 2.1meV!B 
Thus, U/JhS « 4.82 and t/JHS « 0.359. For these pa- 
rameters, the Wigner lattice at low temperatures is stable 
in the phase diagram in Fig. 2. Therefore, the superex- 
change process in Fig. 1(a) should play an important role 
in driving electrons to form the Wigner lattice no matter 
whether the direct nearest neighbor Coulomb interaction 
is strong. It is worth mentioning that the direct Coulonpbi 
interaction will always favor to form the Wigner lattice.ta 
If the direct Coulomb interaction is also included in the 
electronic model, which is not much screened, the stabil- 
ity of Wigner lattice will be greatly enhanced. Note that 
the Coulomb interaction is independent of the magnetic 
structure, and should not be very sensitive to an external 
magnetic field. The effect of field-induced melting of the 
Wigner lattice suggests that the physical origin of the 
state may be closely related to the magnetic structure, 
which is an essential ingredient of the present theory. In 
the actual compounds, both the mechanisms should have 
important impact on the electronic behaviors. It is un- 
likely that only one of them is predominant. As for the 
mean field approximation, when it is sure that the insta- 
bility of Wigner lattice occurs at low temperatures, it is 
an efficient and powerful tool to determine the phase di- 
agram, although some other physical quantities, such as 
critical exponents, cannot be obtained accurately. Due to 
the strongly correlations of electrons, it is still lack of nu- 
merical results to verify the present theoretical prediction 
as this is the first time to discuss instability of the Wigner 
lattice in a model without nearest neighbor or long-range 
interactions. When the system is deviated from x = 0.5, 
the superexchange interaction is still very important to 
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determine the behaviors of electrons. Recently, it was. 
observed that the charge stripes in (La,Ca)Mn03 pair.Ej 
However, the two pairing stripes of Mn'^+ ions are sep- 
arated by a stripe of Mn'*+ ions. This fact suggests the 
nearest neighbor interaction should be very strong. Of 
course, for a comprehensive understanding of the phase 
diagram, including anisotropic properties of charge and 
magnetic orderings, we need to take other effects into 
account. 

The role of Hund's rule coupling in the doped man- 
ganites has been emphasized since the double exchange 
mechanism was proposed. However, the rich phase dia- 
grams in the doped manganites go beyond the picture. 
Our theory shows that the on-site Coulomb interaction 
also has an important impact on the physical properties 
of the system. In the model we investigate, the sign of 
the effective interaction in Eq.(^a|) depends on the ratio 
of jh/U. Repulsive or attractive interaction will leajd to 
quite different physics. In one of our recent papersO, we 
proposed a mechanism of phase separation based on the 
attractive interaction caused by the virtual process (a) in 
Fig. 1, and neglect the on-site interaction U. The phase 
separation can occur in the high and low doping regions. 
As the mechanism of the phase separation is completely 
opposite to the mechanism of the Wigner lattice we dis- 
cuss in this paper, we have to address the issue which 
one occurs for the doped manganites. From the estima- 
tion of the model parameters for the actual compound, 
U / ih ~ 4.82. Thus, the effective interaction should be 
repulsive, not attractive. Frcap this sense, the phase sep- 
aration we predicted in Ref.u could not occur in doped 
manganites. In fact, both the phase separation and the 
Wigner lattice were observed in the family of samples 
with different dopings. For example the phase separatioa. 
was observed in Lai_a;Cua;Mn03 with x=0.05 and 0.08.E2I 
It is worth pointing out that the electronic model is a 
simplified model for doped manganites since the degen- 
eracy of Gg electrons and the Jahn- Teller effect have been 
neglected. The importance of the orbital degeneracy of 
eg electron has been extensively discussed, especially for 
the ferromagnetism near x — Q. If we take into account 
the orbital degeneracy, there may exist an superexchange 
virtual process in the ferromagnetic or A-type antiferro- 
magnetic background, in which the superexchange cou- 
pling between different orbits instead of the spin indices 
in Fig. 1 could produce an attractive interaction as we 
predicted in Ref.l3. The mechanism for phase separation 
may still be responsible for the experimental observation. 
The investigation along this direction is in progress. 

Before ending this paper, we would like to address the 
stability of the Wigner lattice with respect to the transfer 
t. Some experimental analysis suggested that a fslatively 
small t would favor to form the Wigner lattice,cZI which 
seems to be unlikely in contradiction with the phase dia- 
gram in Fig. 2. In the present theory, the Wigner lattice 
occurs in a moderate value of t. On one hand, a large 
t jh), of course, will lead to the instability of the 
Wigner lattice and destroy the double exchange ferro- 



magnetism. In the case, a paramagnetic phase should 
be favored at low temperatures. The perturbation tech- 
nique used in this paper is also not valid. So the region 
of Wigner lattice in Fig. 2 cannot be naively extended to 
the large t case. On the other hand, when t becomes very 
small comparing with jh, the Wigner lattice should also 
be unstable since a small t means to enhance the ratio 
jh/t and a larger ratio is favorable to double exchange 
ferromagnetism. If the antiferromagnetism from t2g elec- 
trons could compete over the double exchange ferromag- 
netism at a; = it would suppress ferromagnetism at 
all the range of xB The effective transfer tcos(0/2) is de- 
termined by either t or the angle of the two spins. The 
Wigner lattice is also accompanied by the strong antifer- 
romagnetic correlation. The field-induced melting effect 
indicates that the Wigner lattice is unstable in the ferro- 
magnetic background, which also indicates the important 
role of the antiferromagnetic correlation to stabilize the 
Wigner lattice. A smaller jaf will reduce the angle 8 and 
should also lead to the instability of the Wigner lattice. 
Thus, a small t does not always favor to form the Wigner 
lattice. 

In short, we derived an effective Hamiltonian for an 
extended Kondo lattice model, based on which a physical 
mechanism for charge ordering in half-doped manganites 
is naturally put forward. 
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